General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



/> IN- FLIGHT ABSOLUTE RADIOMETRIC CALIBRATION OP THE THEMATIC MAPPER 

\ 


K.R. Castle, R.G. Holm, C«7. Kastner, J.M. Palmer, and P.N. Slater 


Optical Sciences Center 
University of Arizona 
Tucson, Arizona 85721, USA 


M. Dinguirard 

Center d'Etudes et de Recherches de Toulouse 
2, Avenue Edouard- Be lin 
31055 Toulouse, Cedex, FRANCE 


C.E. Ezra and R.D. Jackson 
Agricultural Research Service 
U.S. Department of Agriculture 
Phoenix, AZ 85040, USA 


R. Savage 

Atmospheric Sciences Laboratory 
White Sands Missile Range, NM 88002, USA 


(£84-10064) IN-FLIGHT ABSOLUTE R ADIQMETBIC N84- 15633 

CALIBRATION OF THE THEMATIC MAPPER (Arizona 
Univ., Tucson*) 10 p HC AQ2/HF A01 CSCL 14B 

Unclas 
G3/43 00064 



2 


The Thematic Mapper (TM) multi spectral scanner system was placed 
into earth orbit on July 16, 1982 as part of NASA's Landsat-4 payload. 
The entire system was calibrated in an absolute sense at Santa Barbara 
Research Center before launch. 1 * 2 During flight, an internal calibration 
system is used to monitor the calibration of the focal plane.^ In order 
to determine temporal changes of the absolute radiometric calibration of 
the entire system in flight, we initiated a program at White Sands, New 
Mexico, on January 3, 1983, to make spec troradiome trie measurements of 
the ground and the atmosphere simultaneously with TM image collections 
over that area. 4 * By entering the measured values in an atmospheric 
radiative transfer program,^ the radiance levels in four of the spectral 
bands of the TM were determined, band 1: 0.45-0.5 2vm, band 2: 0.53-0.6 lum, 
band 3; 0.62-0.70ym and band 4: 0.'; u-0.91ym. These levels were compared 
to the output digital counts from the detectors that sampled the 
radiometrically measured ground area, thus providing an absolute 
radiometric calibration of the entire TM system utilizing those 
detectors. 

On 3 January 1983, an 80 mm layer of two-day-old snow covered the 
flat gypsum surface at the White Sands Missile Range. The reflectance 
of the snow was measured by reference to a 1.2 m square barium sulfate 
panel using a Barnes Modular Multispectral 8-Channel Radiometer, that 
collected radiant flux simultaneously in all the TM spectral bands over 
a total field angle of 15°. The instrument was mounted on a rotatable 
boom 2.5 m above the snow to allow an average reflectance value to be 
determined for an area of about 0.5 x 5 m. The measurements were made 
at 10.08 hours coinciding with the overpass of the TM- and a solar zenith 
angle of 60°. When ratioed to the reflectance of the barium sulfate 
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panel at 60*, the reflectance of the snow in TM bands 1 to A was found 
to be 0.769, 0.761, 0.756 and 0.732 respectively, with an Yas uncertainty 
of £0.023 in all cases. 

A solar radiometer^* using nine 10-nm spectral bands in the visible 
and near ir was used to determine the total spectral optical depths t t 
at these nine wavelengths. The measured barometric pressure of 689.5 mB 
allowed the Rayleigh spectral optical depth T a ,y to be determined. From 
these data the Mie and ozone spectral optical depths were determined. 

In bands 1, 1, and 3, the component of due to molecular 
absorption. Tabs, is entirely due to ozone. Water vapor and CO 2 are 
present in. addition to ozone in band 4. Their effects have been included 
in Table 1 which lists the values of the various atmospheric 
constituents in bands 1 through 4. 


Table 1. Atmospheric values. 


Band 

* c (um) 

t t 

T Mie 

* T Ray 

T abs 

1 

0.485 

0.2913 

0.1475 

0.1424 

0.0014 

2 

0.57 

0.2177 

0.1382 

0.0736 

0.0059 

3 

0.66 

0.1 /i.6 

0.1282 

0.0406 

0.0030 

4 

0.84 

0.1344 

0.1099 

0.0153 

0.0092 

Using 

the data 

in Table 

1 as input 

to the radiative 

transfer 


following quantities were calculated: 

Eq Dir : the direct solar irradiance at the ground * cosG z exp(-x^ sec@ z ) 

Ed DifS the diffuse solar irradiance at the ground 

L Dir s t * ie direc t radiance at the TM due to J^Dir + ^D,Dif 

Lp: the path radiance at the TM * Lj - (Ej^Dir + E D,Dif^ f - exp(-x^ sec5 # ) 

Lf! the total radiance at the TM 
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Table 2. I r radiance 3 and radiances (normalized to unity solar 
exoatmospheric irradiance). 


Solar 

Zenith 


Band 

Angle 

E D,Dir 

E D,Dif 

L Dir 

L P 

4 

1 

55* 

0.345 

0.185 

0.097 

0.033 

0.1301 


65* 

0.212 

0.152 

0.067 

0.025 

0.092 

2 

55* 

0.392 

0.145 

0.105 

0.024 

0.129 


65* 

0.253 

0.122 

0.043 

0.019 

0.092 

3 

55° 

0.425 

0.122 

0.111 

0.019 

0.130 


65* 

0.282 

0.104 

0.078 

0.015 

0.093 

4 

55* 

0.454 

0.095 

0.112 

0.014 

0.126 


65* 

0.308 

0.082 

0.079 

0.011 

0.090 


Using the equivalent passband technique of Palmer and Tomasko^ 
and the exoatmospheric solar spectral irradiance values of Neckel and 
Labs,** the values for the exoatmospheric irradiances within the TM 
passbands (ExAA), were calculated for the earth-sun distance on 
January 3. The required values of Lj. for the solar zenith angle of 
60° were found by interpolating between the 55° and 65° data from 
Table 2. These, when multiplied by their corresponding Ex&A values, 
gave the radiances in mW/cm^~sr in the TM passbands listed in Table 3. 
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Table 3. Exoatmospheric irradiances and the radianeea at the TM in 
the TM paaabands. 


Band 

Equivalent 
bandwidth in pm 

ExBX 

mW/cm^ 

W 2 

mW/cm -sr 

1 

0.0712 

12.2 

1.36 

2 

0.0884 

13.8 

1.52 

3 

0.0773 

10.7 

1.19 

4 

0.1345 

13.7 

1.47 


By identifying our site on the raw image data, we determined 
which detectors scanned the area and how many samples each collected. 
Detectors l, 2, 3, 15 and 16 collected 5, 3, l, 2, and 4 samples 
respectively. We found the average digital count for each detector, 
then calculated the average detector spectral radiance using offset 
and gain values reported by Barker et al.* These spectral radiances 
were multiplied by the number of samples for each detector. The 
resultant products were added and then divided by the total number of 
samples, 15. Thus we derived a value for the average spectral 
radiance of our site as measured by the TM, proportionally weighted 
according to the number of samples per detector. 

The fina 1 step in the calculation was to multiply the average 
spectral radiance at the TM by the equivalent passband to determine 
the radiance within each TM band (Table 3). These values are compared 
in Tables 4, 5 and 6. 
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Table 4. Calculation of radiance In TM band 2 from preflight 
calibration data 1 * and comparison with inflight White Sands data. 


Detector 

Average 

Average Spectral 

Count Radiance 

Gain 

Offset 

3 

139.0 

17.04 

8.014 

2.41 

2 

143.7 

17.42 

8.117 

2.33 

1 

146.2 

17.52 

8.174 

2.99 

16 

141.5 

17.43 

7.979 

2.43 

15 

147.5 

17.71 

8.195 

2.83 

Weighted average 
Equivalent band 

- 17.47 
radiance 

mW cm”2 ar -l 

from weighted average “ 1.537 

mW cm * sr’ 


Agreement » -1.0 % and +2.5 % with reference to preflight calibration 
and December 8, 1982 inflight data^ respectively. 


Table 5. Calculation of radiance in TM band 3 from preflight 
calibration data* and comparison with inflight White Sands data. 


Detector 

Average 

Count 

Average 

Spectral 

Radiance 

Gain 

Offset 

3 

165.0 

15.40 

10.590 

1.89 

2 

169.7 

15.86 

10.602 

1.57 

1 

171.6 

15.73 

10.777 

2.13 

16 

167.8 

15.85 

10.484 

1.63 

15 

172.5 

15.88 

10.769 

1.53 

Weighted average 
Equivalent band 

■ 15.78 mW cm 
radiance from 

“2 s r "l um'^ 
weighted average • 1.215 

mW cm" sr' J 

t 


Agreement ■ -1.8% and 3.0% with reference to preflight and December 
8, 1982 inflight data^ respectively. 
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Table 6. Calculation cf radiance In TM band 4 from preflight 
calibration data*- and comparison with inflight White Sands data. 


Detector 

Average 

Count 

Average 

Spectral 

Padiance 

Gain 

Offset 

3 

133.0 

11.89 

11.019 

1.96 

2 

132.0 

12.02 

10.817 

1.94 

1 

134.2 

12.00 

10.972 

• 2.53 

16 

132.0 

12.01 

10.828 

2.00 

15 

131.5 

12.05 

10.771 

1.69 


Weighted average • 12.01 mWcnT^ sr’ 1 Pm’ 1 
Equivalent band radiance from weighted average * 
Equivalent band radiance from Table 3 * 1.47 mWcm 


1.609 mWcm* 
’ 2 sr’ 1 


sr 


-1 


Agreement * -8.5% and -5.6% with reference to preflight and December 
8, 1982 inflight data 2 respectively. 


Because of the litiitations of the instruments used on January 3, 
T a j )S in band 4 due to water vapor and CO 2 could not be measured. 
Instead, the LOWTRAN code was used to compute the water vapor and CO 2 
transmittances across band 4 at 5 cm* 1 intervals. The water vapor 
transmittance was scaled, to account for the 44% relative humidity 
measured at White Sands on January 3, then averaged to find 
Finally, the measured value of t oz * 0.0003 was added to provide T a k s . 
The effect of an error in this estimate can be judged by noting that 
the inclusion of water vapor and CO 2 lowered the predicted radiance 
level at the TM by only 2%. 

Band 1 saturated over the snow field at White Sands. Pre-flight 
data 1 indicate that a saturation level of 255 counts corresponds to a 
radiance at the sensor of 1.14 mWcm'^sr’ 1 in TM band 1. Our estimate 
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Is that the snowfield provided a radiance level of 1.36 mUcm~^ar~^ at 
the sensor. 

The estimated uncertainty In the results of this first 
measurement is ±5%, while the estimated uncertainty in the preflight 
calibration is no better than ±6%.^*® Ve are encouraged that our 
results for TM bands 2, 3 and 4 fall between 1% and 8.5% of the 
preflight measurements and between 2.5% and 5.6% of the December 8, 
1982 inflight data. 

Ue are presently fabricating field equipment to provide more 
detailed and accurate measurements of the surface and atmosphere at 
Uhite Sands. Our goal is to reduce the uncertainty in sensor absolute 
calibration to less than ±3%. Ue plan to continue the work described 
here to include the inflight absolute radiometric calibration of the 
second TM and the Systfcme Probatoire d'Observation de la Terre, Haute 
Resolution Visible (SPOT/HRV) systems. 

We wish to thank B.M. Herman for the use of his radiative 
transfer code and J.A. Reagan for the use of his solar radiometer. The 
University of Arizona authors wish to acknowledge support from NASA 
contract number NAS5-27382, J.L. Barker, Science Representative, and 
the USDA, broadform contract number 12-14-5001-38. 
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